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We investigate the hole spin relaxation due to the Rashba spin-orbit coupling induced by an 
external perpendicular electric field in bilayer WSe 2 . The Rashba spin-orbit coupling coefficients 
in bilayer WSe 2 are constructed from the corresponding monolayer ones. In contrast to monolayer 
WSe 2 , the out-of-plane component of the bilayer Rashba spin-orbit coupling acts as a Zeeman- 
like held with opposite directions but identical values in the two valleys. For in-plane spins, this 
Zeeman-like held, together with the intervalley hole-phonon scattering, opens an intervalley spin 
relaxation channel, which is found to dominate the in-plane spin relaxation in bilayer WSe 2 even 
at low temperature. For out-of-plane spins, this Zeeman-like held is superimposed by the identical 
Hartree-Fock ehective magnetic helds in the two valleys, and hence diherent total effective magnetic 
helds between two valleys are obtained. Owing to the large difference of the total helds at large spin 
polarization, diherent out-of-plane spin relaxation times in the two valleys are obtained when the 
intervalley hole-phonon scattering is weak at low temperature and low hole density. This diherence 
in the spin relaxation times can be suppressed by enhancing the intervalley hole-phonon scattering 
through increasing temperature or hole density. Moreover, at large spin polarization and low tem¬ 
perature, due to the weak intravalley hole-phonon scattering but relatively strong hole-hole Coulomb 
scattering, the fast spin precessions are found to result in a quasi hot-hole Fermi distribution char¬ 
acterized by an ehective hot-hole temperature larger than the temperature, which also enhances the 
intervalley scattering. During this process, it is interesting to discover that the initially equal hole 
densities in the two valleys are broken in the temporal evolution, and a valley polarization is built 
up. It is further revealed that this comes from the diherent spin relaxation processes at large spin 
polarization and diherent spin-conserving intervalley scattering rates between spin-up and -down 
holes due to the diherent ehective hot-hole temperatures. 

PACS numbers: 72.25.Rb, 71.10.W, 71.70.Ej, 72.10.Di 


I. INTRODUCTION 

Very recently, the spin dynamics in monolayer (ML) 
transition metal dichalcogenides (TMDs) has attracted 
much attention,i“— partly due to their unique electrio^^— 
and novel optica h^d^~ — properties. Specifically, owing 
to the direct gap at the K (K') pointer— and large en¬ 
ergy spin splitting of the valence bands /^d^d’^ the chi¬ 
ral optical valley selection rule in ML TMDsi^Ti^^ al¬ 
lows optical control of the valley pseudospin^^^— and real 
spinfiip— making them promising candidates for the 
spintronic application. 

Among these studies in ML TMDs, carrier spin relax¬ 
ation due to the Elliot- Yafe t^^d^ (EY) and D’yakonov- 
Perel— (DP) mechanisms is an important property to 
understand for any possible spintronic applications. For 
the intrinsic EY mechanism, it was claimed that the in¬ 
tervalley out-of-plane spin-flip scattering is forbidden by 
the time reversal symmetry, while the flexural phonon 
vibrations thermally activated at high temperature can 
lead to the intravalley out-of-plane spin relaxation.-!- Ad¬ 
ditionally, due to the marginal in-plane spin mixing, the 
contribution of the EY mechanism to the in-plane spin 
relaxation is negligible^ Eor the DP mechanism, in the 
intrinsic situation, the out-of-plane spin relaxation pro¬ 
cesses are forbidden due the good quantum number of 
az but the in-plane processes exist including the in¬ 


travalley and intervalley processes. Particularly, the in¬ 
tervalley process arises from the intrinsic spin-orbit cou¬ 
pling (SOC), which acts as opposite effective magnetic 
fields in the two valleys and hence opens an intervalley 
in-plane spin relaxation channel in the presence of the 
intervalley carrier-phonon scatterings. The extrinsic in¬ 
fluence such as the flexural deformations;^ the external 
in-plane magnetic field^ or the Rashba SOC induced by 
an external perpendicular electric field;^ can cause the 
out-of-plane spin relaxation. Moreover, in contrast to 
the electron spin relaxation, the hole spin relaxation pro¬ 
cesses are markedly suppressed due to the large energy 
splitting of the valence bands 

Consisting of two-layer TMDs, bilayer (BL) TMDs 
also obey the chiral optical valley selection rule and 
possess good spin charactersapart from new fea¬ 
tures such as layer pseudospin,-^ electrical-tuned mag¬ 
netic moments^ and magnetoelectric effect)^ Specifi¬ 
cally, due to the 180° in-plane rotation between the up¬ 
per and lower layers in BL TMDs, in the K (K') valley, 
the hole bands in one layer have opposite spin polar¬ 
izations between the energy-degenerate ones in another 
layer. The interlayer hopping, which exists only in holes 
in the same valley with the same spin, has no influence on 
this spin degeneracy. This is very different from the ML 
TMDs, where the hole bands possessing opposite spin po¬ 
larizations in the K (K') valley are largely energy split. 
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markedly suppressing the hole spin relaxation. In BL 
TMDs, with the intrinsic EY hole spin relaxation and 
the intrinsic DP in-plane one suppressed in each layer 
while the intrinsic DP out-of-plane one forbidden, the 
Rashba SOC of the lowest two degenerate hole bands^ 

rj'" = [ - v{l + ak‘^)ky,i^{l + ak^)kx, fJ.r]]Ez, (1) 

induced by an external perpendicular electric field E^, 
is expected to make the dominant contribution to the 
hole spin relaxation. Here, the z-axis is set to be per¬ 
pendicular to the BL TMD plane; ^ = 1(—1) represent¬ 
ing the K (K') valley. The out-of-plane component of 
Rashba SOC serves as a Zeeman-like term nijEz with 
opposite effective magnetic fields in the two valleys j^^i^^ 
Due to this coupling of real spin and valley pseudospin, 
one may expect the interplay of the spin polarization 
with the valley polarization and the rich physics of spin 
and valley dynamics in BL TMDs. Specifically, for in¬ 
plane spins, this Zeeman-like term, together with the 
intervalley hole-phonon scattering, opens a new inter¬ 
valley spin relaxation channel. For out-of-plane spins, 
this Zeeman-like term is superimposed by the Coulomb 
Hartree-Fock (HF) self-energy, which has been under¬ 
stood first theoretically^^— and then experimentally^ 
in semiconductors to serve as an effective magnetic field 

fiHF(k) =-^Vk-k'Tr[pk^cr], (2) 

with Vk-k' being the screened Coulomb potential. In the 
case of the valley-independent out-of-plane spin polariza¬ 
tion, the total effective magnetic fields 

= (^hf + mEz)ez (3) 

have different values in the two valleys, leading to the 
valley-dependence of the out-of-plane spin relaxation. 

In the present work, by utilizing the kinetic spin 
Bloch equation (KSBE) approach^ with the hole-hole 
Coulomb, (both the intra- and intervalley) hole-phonon, 
and long-range hole-impurity scatterings included, we in¬ 
vestigate the hole spin relaxation in BL WSe 2 due to 
the Rashba SOC [Eq. ([T])]. In our investigation, the ini¬ 
tial occupations of holes of each spin are identical in the 
two valleys. In the case of small spin polarization (i.e., 
weak HF effective magnetic field), the total effective mag¬ 
netic fields [Eq. ((S])] determined by the Zeeman-like fields, 
have identical absolute values in the two valleys. We in¬ 
vestigate the temperature and hole density dependence 
of both the out-of- and in-plane spin relaxation times 
(SRTs), and show that for in-plane spins, the relaxation 
process is dominated by the intervalley spin relaxation 
channel induced by the Zeeman-like fields, and the contri¬ 
bution of the intervalley hole-phonon scattering becomes 
dominant. But for out-of-plane spin relaxation, the in¬ 
tervalley hole-phonon scattering is marginal, which indi¬ 
cates that out-of-plane spins relax independently in the 
two valleys. With the intervalley scattering removed, the 


out-of-plane SRTs in the strong scattering limi t^^'^^ 

- (Di(k))r, (Hi(k))r, ^ (Di(k)) ^ ^ 

are identical in the two valleys due to the same absolute 
values of the total effective magnetic fields. Here, r2j^(k) 
and Tp are the in-plane Rashba terms and momentum 
relaxation time due to the intravalley scattering, respec¬ 
tively. 

In the case of large spin polarization (i.e., strong HF 
effective magnetic field), for in-plane spins, the SRT is 
insensitive to the spin polarization, similar to the case 
of ML MoS 2 ^ For out-of-plane spins, a large difference 
of the total effective magnetic fields [Eq. ([3])] between 
the two valleys is obtained. We find that the intervalley 
hole-phonon scattering in this situation is weak at low 
temperature and low hole density, and hence the out-of- 
plane spins relax independently in the two valleys, lead¬ 
ing to different SRTs [Eq. (|4])]. The enhancement of the 
intervalley hole-phonon scattering by increasing temper¬ 
ature can suppress this difference in the SRTs. Moreover, 
the intervalley scattering also becomes stronger with the 
increase of hole density at low temperature, which origi¬ 
nates from the different SRTs. Specifically, as shown in 
Fig.jTJ when the HF effective magnetic field and Zeeman- 
like field have opposite directions, say, in the K valley, the 
smaller total effective magnetic field [Eq. dS])] leads to a 
faster SRT [Eq. (jj])] in this valley. Over time, this faster 
spin relaxation makes the density for spin-down (-up) 
holes larger (smaller) than the corresponding one with 
the same spin in the K' valley, triggering the intervalley 
scattering of spin-down (-up) holes from the K (K') val¬ 
ley to the K' (K) one by emitting phonons to suppress 
this density difference. With larger density difference by 
increasing the hole density, the intervalley scattering be¬ 
comes stronger. In addition, due to the weak intravalley 
hole-phonon scattering at low temperature but relatively 
strong hole-hole Coulomb scattering, we find that the fast 
spin precessions result in a quasi hot-hole Fermi distri¬ 
bution characterized by an effective hot-hole temperature 
Teff larger than the temperature T, which also enhances 
the intervalley scattering. 

During above process, it is interesting to discover that 
the initially equal hole densities in the two valleys are 
broken in the temporal evolution, with more holes are 
accumulated in the K' valley, leading to the build up of 
the valley polarization. This arises from the larger effec¬ 
tive hot-hole temperature for spin-down holes than that 
for spin-up ones, which makes the spin-conserving inter¬ 
valley scattering rate of spin-down holes faster 

than that of spin-up holes (see Fig. [T]). Specif¬ 

ically, due to the large difference of Fermi energies of 
spin-up and -down holes, the Rashba SOC [Eq. (P)] near 
the Fermi energy of spin-up holes transfers holes from 
the spin-up states into the spin-down ones with the same 
energies, and hence more spin-down holes occupying the 
states with the energies higher than the corresponding 
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Fermi energy, leading to the larger effective hot-hole tem¬ 
perature. According to our calculation, with the exper¬ 
imental obtainable hole density and the spin polariza¬ 
tion reaching 60 %, the accessible valley polarization can 
reach beyond 1 % which can last hundreds of picosec¬ 
onds. In addition, due to the absence of electrons in 
the conduction bands in p-type BL WSe 2 , the valley- 
depolarization induced by the exchange interaction^ is 
absent here. Therefore, the valley polarization proposed 
above can be measured experimentally. This has not yet 
been reported in the literature. 


K K' 

It tt 



FIG. 1: (Color online) Schematic of the different spin relax¬ 
ation processes in the two valleys and valley polarization pro¬ 
cess. In the figure, | 'll (fl) uh(ih)} stands for the upper- (lower-) 
layer states for spin-up (-down) holes; the purple (gray) filled 
arrows, which have the same (opposite) directions in the two 
valleys, denote the HF (Zeeman-like) effective magnetic fields; 
the dashed lines stand for the initial occupations of holes of 
each spin. On one hand, this schematic shows that due to the 
smaller total effective magnetic field in the K valley, the SRT 
is faster than , which makes the density for spin-down 
(-up) holes in the K valley larger (smaller) than that in the 
K' valley in the temporal evolution. On other hand, the in¬ 
tervalley scattering time of spin-down holes is faster 

than that of spin-up holes due to the larger effective 

hot-hole temperature for spin-down holes (red). 

This paper is organized as follows. In Sec. m we in¬ 
troduce our model and lay out the KSBEs. Then in 
Sec. IIII A1 we investigate the temperature and the hole 
density dependence of the out-of-plane spin relaxation in 
BL WSe 2 with weak HF effective magnetic field. The 
anisotropy of the spin relaxation is also presented in this 
part. In Sec. IIII Bl we show the influence of HF effective 
magnetic field when the large out-of-plane spin polariza¬ 
tion is considered. The valley polarization due to the 
large spin polarization is also addressed in this part. We 
summarize in Sec. El 


II. MODEL AND KSBES 

In BL TMDs, in the presence of an external perpendic¬ 
ular electric field, the effective Hamiltonian of the lowest 


two hole bands near the K (K') point can be written as^ 

H^s = • < t / 2 , ( 5 ) 

where e^k = k^/(2m*), with m* being the effective mass 
of the hole; the Rashba SOC is given in Eq. ©■ How¬ 
ever, the Rashba coefficients zz, a and r] are absent in the 
literature. In this work, we use the ML Rashba SOG^ in 
each layer to construct the BL Rashba SOC. Specifically, 
a hole in the lowest hole band in a given layer can first 
hop into the second layer by the spin-conserving inter¬ 
layer hopping and then flip spin by the ML Rashba SOC 
to the lowest hole band in the second layer. Or the hole 
in the lowest hole band can first flip spin in the first layer 
and then hop to the second layer with same spin. The 
lowest two hole bands in BL TMDs, which have oppo¬ 
site spin polarizations, are then coupled with each other. 
Due to the large spin splitting of the hole bands in each 
layer, both the ML Rashba SOC and interlayer hopping 
can be treated perturbatively. We use the Lowdin par¬ 
tition method^i^ to derive the effective Rashba SOC of 
the lowest two hole bands in BL TMDs, and obtain the 
Rashba coefficients, which are given in Appendix in 
details. 

The KSBEs are then written 

dtPfilc — dtPfiU.\cQh ff dtPfiU.\ scat j (6) 

where represent the density matrices of holes with the 
off-diagonal terms p^k.i-i = Pj^k-ii representing the 

spin coherence and the diagonal ones Pfj.k.a-(T=fiJ.k,tT {o' = 
±^) being the hole distribution functions. The coherent 
terms 9tp^k|coh describe the spin precessions of holes due 
to the effective magnetic field and the HF effective 
magnetic field Hhf- Tbe scattering terms dtp^ulscat in¬ 
clude the hole-hole Coulomb, long-range hole-impurity, 
intravalley hole-phonon and especially the intervalley 
hole-phonon scattering. The detailed expressions for the 
coherent and scattering terms can be found in Ref. Is^ . 
Particularly, the large Zeeman-like term in the 

Rashba SOC [Eq. ([T])], leads to the energy-splitting of 
the lowest two hole bands in BL TMDs, which has oppo¬ 
site signs in the two valleys, and hence affects the spin- 
conserving intervalley hole-phonon scattering but has no 
effect on the intravalley scattering. By considering this 
energy-splitting, the intervalley scattering parts in the 
KSBEs are given by^ 

^zPMkl^r = >) + 5^k(<, >)t} - {<0>}, (7) 

with 

Sfik{>,<)=TT ^ |-^/i/i'qPP^'k-q^iLl''7l^M'Z2P[rk[-^q 

qMOl'72 

xd(£jy2k ^r^ik— q ff ^q) T -^q ^(^T/ik— q ^?72k ^q)] ■ (^) 

Here, = Aq is the phonon distribution, = Aq -|- 

PjTk = PMk> and p>^ = I - p^k; Piim) = ±1 and 
e±ik = Ek ± pEz] The projector operator reads = 
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(1 + ?7^(Tz)/ 2; is the scattering matrix element 

of the intervalley phonon mode A. 

However, the hole-phonon scattering matrix elements 
in BL TMDs have not yet been reported in the literature. 
For the intervalley hole-phonon scattering, with the inter¬ 
valley scattering in each layer suppressed due to the large 
spin splitting holes in a given valley at a given layer are 
scattered into the other valley at different layer. We de¬ 
rive the matrix elements of this intervalley scattering by 
using the tight-binding model according to the arXiv ver¬ 
sion of the work by Viljas and Heikkilai^ We show that 
only the phonon modes and make the main con¬ 
tribution. Here, Kq {K^) is the phonon mode at the 
K point corresponding to the irreducible representation 
E2 of group Csh with the lower (higher) phonon energy.— 
For the intravalley hole-phonon scattering, we derive the 
matrix elements of out-of-plane phonon according to the 
same work by Viljas and Heikkila; the matrix elements of 
in-plane phonon in BL TMDs are constructed by using 
the ML ones, which have been reported in the work by 
Jin et al.^ The specific derivation of the hole-phonon 
scattering matrix elements is give in Appendix [C] 

The matrix elements of the hole-phonon scatter- 

rL 

ings including the intervalley K\ and 

rH 

phonon scatterings and the intravalley in¬ 
plane acoustic (AC) (|Af^qp), in-plane optical (OP) 

mlXW and out-of-plane OP 

gi 

phonon scatterings are given by 


|2 _ 4 


'"'“-‘‘■I 

l<;,t = 


I flfl'Ci 


(S )4 , 

2 pVLA 


9 ,/ 2 
.B2g |2 _ 4 j 




piq 


,2 




P^T,E^^ 


42 d1 

|M 

I /i/i'q I 


t 


/ 2 




2M, 


Mt 




( 9 ) 

( 10 ) 

( 11 ) 

( 12 ) 

( 13 ) 

( 14 ) 


where A|g, are the phonon 

modes at the F point in BL TMDsj^^— Hr(K),A is the 
energy of phonon mode A at F(K) point; p is the mass 
density; v-^a is the sound velocity corresp onding to lo n- 
gitudinal acoustic (LA) phonon; 7 = ^JUbT/{ pv^) is 
the parameter for the energy dispersion of the out-of- 
plane AC phonon with vq being the corresponding sound 
velocityji ks represents the Boltzmann constant; t'j_ is 
the derivative of interlayer hopping t± with respect to 
the corresponding hopping bond length; S and Dop are 
the deformation potentials of the in-plane AC and OP 
phonon in ML TMDs;^ respectively; Mm and Mx are 


the masses of the M atom and the X atom of the TMD 
MX2, respectively. The remaining matrix elements of the 
hole-hole Coulomb and long-range hole-impurity scatter¬ 
ings can be found in the work by Wang and Wu.- 


III. NUMERICAL RESULTS 


TABLE I: Parameters for WSe 2 used in our calculation. Note 


that mo stands for the free electron mass. 


uiM) 

0.0342“ 

Tj (A) 

1.06“ 

m* j mo 

0.51® 

a(A)^ 

-5.673“ 

K 

5.20“ 

4 (eV/A) 

1.26"^ 

H(eV) 

2.1“ 

Dop (10® eV/cm) 3.1“ 

vla (10® cm/s) 

3.30“ 

I’D (10^ cm/s) 

2.24“ 

(meV) 

30.8“ 

(meV) 

38.51-^ 

(meV) 

17.5“ 

^k,kh (meV) 

31.5“ 

Mse/Mw 

0.429 

P (10"' g/cm^) 

6.6“ 

“ Appendix 

A. ® RefJ^. 

“ Ref. [la 

“ Re^. 

^ Ref. Ha 

Ref. 113 . 


In this work, we focused on the case of WSe2. The ini¬ 
tial occupations of holes of each spin are identical in the 
two valleys in our calculation. The spin-polarization di¬ 
rection is along the z-axis unless otherwise specified. The 
long-range impurity density is taken to be W = 0 . 02 V;, 
with Nh being the hole density. Then the correspond¬ 
ing mobility at T = 250 K in our investigation is around 
180 cm^/(V's), which agrees with those reported in the 
existing experiments<^“— All the material parameters 
used in our calculation are listed in Table H) With above 
parameters, by numerically solving the KSBEs, we dis¬ 
cuss the spin relaxation at small and large spin polariza¬ 
tions (i.e., weak and strong HF effective magnetic fields) 
in Secs. IHI K\ and HIT^ respectively. 


A. Weak HF effective magnetic field 

In this subsection, the initial spin polarizations Pg in 
the two valleys are set to be 2.5 %, leading to the weak 
HF effective magnetic field (|Hhf|"C|77Aj;|) in our calcula¬ 
tion. Then the out-of-plane SRTs due to the intraval¬ 
ley scatterings [given in Eq. (| 1 ]) with Hgg Ri ppEz], have 
identical values Tgz in the two valleys. Moreover, ow¬ 
ing to the competition between [|Hggp/(H5_(k))]Tp and 
[(r 25 _(k))rp]“^ in Eq. (g]), the SRT can be divided into two 
regimes: regime I, the anomalous EY-like regime, where 
the first term makes a more important contribution and 
hence the SRT shows the EY-like behavior; regime H, 
the normal strong scattering regime, where the second 
term becomes more important. The crossover between 
regimes I and H is determined by |Heffrp|=l with the 
system sitting in regime I (H) when jOeffTpl > (<) 1 . In 
this work, with jHeffj ~ \vEz\ at small spin polariza¬ 
tion, we discuss the temperature and hole density depen¬ 
dence of the out-of-plane SRT at the small {\riEzTp\<^l) 
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and medium {\'qE^Tp\~l) fields. The large field with 
\riEzTp\:>l, which leads to a large energy-splitting of the 
hole bands similar to ML WSe 2 , is not included in this 
work. We also investigate the anisotropic spin relaxation 
by varying the spin-polarization direction in this subsec¬ 
tion. 


1. Temperature dependence of out-of-plane spin relaxation 

We first investigate the temperature dependence of the 
out-of-plane spin relaxation. The SRTs Tgz as function of 
temperature T at different hole densities are plotted in 
Figs. [5] and [3] at the small [E^ = 0.003 V/A) and medium 
{Ez = 0.01 V/A) helds, respectively. All the relevant 
scatterings are included in the calculation. From Fig. [21 
it is noted that the intervalley hole-phonon and the long- 
range hole-impurity scatterings are marginal to the out- 
of-plane spin relaxation, since the SRTs with both the in¬ 
tervalley and long-range hole-impurity scatterings (curve 
with open circles) removed are close to the one with all 
scatterings included (curve with squares). To further 
elucidate the role of the remaining hole-hole Coulomb 
and intravalley hole-phonon scatterings, we compare the 
SRT with only the hole-hole Coulomb scattering included 
(curve with crosses) and the one with only the intraval¬ 
ley hole-phonon scattering included (curve with dots). 
The hole-hole Coulomb scattering is found to be more 
important at low temperature since the SRT with only 
the hole-hole Coulomb scattering is much larger than the 
counterpart, while the intravalley hole-phonon scattering 
plays an important role at high temperature. 

At the small field (the SRT is in regime II), a 
peak around T = 60 K at the hole density Nh = 
0.4x10^^ cm“^ (the Fermi temperature Tp = 11 K) is 
observed in the temperature dependence of the SRT [see 
Fig. |2ta)]. As revealed in the previous works j^^'^^i^^ 
this peak is due to the crossover of the hole-hole 
Coulomb scattering from the degenerate to nondegen¬ 
erate limits [i.e., l/Tp'^ix.h\{Tp/T)T"^/Tp when T<^Tp 
and l/rp''ocl/T when T^Tp] ^^^™ when the SRT is 
in the normal strong scattering regime. Unlike the 
previous studies in semiconductors, where the peak 
locations in different materials locate in the range 
{Tp/4:,2Tp)^^^L^^-M the peak location of p-type 
WSe 2 is around 5Tp. Moreover, compared with the 
case with only the hole-hole Coulomb scattering, with 
all scatterings included, the peak is shifted toward a 
higher temperature at for the low density and found 
to be destroyed (not shown) when 7V;i>2xl0^^ cm“^ 
(5 Tf> 275 K). This is due to the fact that the reduction 
of the Coulomb scattering in the nondegenerate limit can 
be compensated by the increase of the intravalley hole- 
phonon scattering, similar to the previous work)^ and 
the intravally hole-phonon scattering is stronger than 
the hole-hole Coulomb scattering after 275 K. Conse¬ 
quently, when 7V;i>2xlO^^ cm“^, the SRT increases with 
the decrease of Tp by increasing temperature, as shown 



FIG. 2: (Color online) The SRTs Tsz as function of T at 
(a) Nh = 0.4x10^^ cm“^, (b) Nh = 0.8x10^^ cm”^, and 
(c) Nh = 2x10^® cm~^. Squares: all the relevant scat¬ 
terings are included; Triangles: the intervalley hole-phonon 
scattering is removed; Open circles: both the intervalley hole- 
phonon and long-range hole-impurity scatterings are removed; 
Crosses (Dots): only the hole-hole Coulomb (intravalley hole- 
phonon) scattering is included. Ez = 0.003 V/A. 


in Fig. |2Kc) at Nh = 2x10^^ cm 

At the medium held (see Fig. [3]), the crossover be¬ 
tween regimes I and II is observed at Tf. Specihcally, 
with the increase of Tp by removing each scattering at 
the same temperature, Tsz increases in regime I while 
decreases in regime II. Compared with the location of 
at the small held, which is determined by Tp and 
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FIG. 3: (Color online) The SRTs Taz as function of T at (a) 
Nh = 0.4x10^^ cni~^, (b) Nh = 0.8x10^^ cni“^, and (c) 
Nh = 2x10^® cni“^. Squares: all the relevant scatterings are 
included; Triangles: the intervalley hole-phonon scattering is 
removed; Open circles: both the intervalley hole-phonon and 
long-range hole-impurity scatterings are removed; Crosses: 
only the hole-hole Coulomb scattering is included. The in¬ 
set in (a) zooms the temperature range 30-31 K. The vertical 
black dashed and red dash-dot lines indicate TJ and T^, re¬ 
spectively. Ez = 0.01 V/A. 


not influenced by E^, we find that the location of the 
crossover between regimes I and II TJ (indicated by the 
vertical black dashed line) is always smaller than (in¬ 
dicated by the vertical red dash-dot line) at the corre¬ 
sponding hole density. This is because that with in¬ 
creasing temperature before the decrease of Tp by 


the hole-hole Coulomb and intravalley hole-phonon scat¬ 
terings leads to the crossover from regime I {\r]EzTp\ > 1) 
to II {IrjEzTpl < I). Therefore, with increasing tempera¬ 
ture before Tj}, the SRT in regime I {tsz oc Tp) decreases 
whereas in regime II {tsz oc t“^) increases [see Fig.[3l[b)]. 
Additionally, with increasing temperature after (the 
SRT is entirely in regime II), as mentioned above, Tp 
increases (decreases) when Nh< (>)2xl0^^ cm“^, and 
hence Tgz decreases (increases), as shown in the fig¬ 
ures. Consequently, when is close to T^, as shown 
in Fig. [3l[a) at Nh = 0.4x10^^ cm“^, the peak existed 
at the previous small field with the same hole density is 
vanished here. But when TJ is separated from T^, as 
shown in Fig. [3][b) at Nh = 0.8x10^^ cm“^, the peak 
existed at the small field with the identical hole density 
[Fig. mb)] is observed at the same temperature position 
here, but is less visible. 

2. Hole density dependence of out-of-plane spin relaxation 



FIG. 4: (Color online) The SRTs Taz versus Nh at differ¬ 
ent temperatures with (a) Ez = 0.003 V/A and (b) Ez = 
0.01 V/A. The vertical red (blue) dashed line indicates the 
crossover between regimes I and II at T = 30 (60) K. 

Next we turn to study the hole density dependence 
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of out-of-plane spin relaxation. The SRTs Tsz versus 
the hole density at different temperatures are plot¬ 
ted in Figs.mja) and (b) at the small {E^ = 0.003 V/A) 
and medium = 0.01 V/A) helds, respectively. At 
the small field, a peak around Tp^T/A.B is observed in 
the density dependence of the SRT, which is due to the 
crossover of holes from the nondegenerate to degener¬ 
ate limits when the SRT is in the normal strong scatter¬ 
ing regimeAt the medium field, similar peak is 
also observed, as shown in Fig. S^b), indicating the SRT 
around TpKiT/A.b still sits in regime IL But compared 
with the results at the small field, it is noted that the de¬ 
crease of the SRT in the degenerate limit at the medium 
field slows down after = 0.8xl0^^cm“^ at T = 30 K 
(curve with squares). This arises from the crossover be¬ 
tween regimes I and II. Specifically, with increasing the 
hole density in the degenerate limit at low temperature, 
the increase of \r]EzTp\ due to the dominant hole-hole 
Coulomb scattering [i.e., I/Tp*'ocln(rF/r)T^/Tp when 
Tp^T {Tp(xNh for the two-dimensional case)] drives the 
SRT from regime II {\riEzTp\ < I) to regime I {\riEzTp\ > 
I). In the degenerate limit, one has TszOC ln{Tp/T)T^/Tp 
in regime II whereas rszOcI/[T^ ln(Tp/T)] in regime I 
since {^l\{'k))(xTp in Eq. (jd]). Consequently, with the 
increase of Nh, the decrease of SRT slows down after 
the crossover. Similar behavior can also be observed at 
T = 60 K (curve with triangles). But at high temper¬ 
ature T = 150 K (curve with crosses), the intravalley 
hole-phonon scattering becomes stronger and the reduc¬ 
tion of the Coulomb scattering in the degenerate limit 
can be compensated by the increase of the hole-phonon 
scattering, leading to the crossover indistinguishable. 


3. Anisotropic spin relaxation 

We also address the anisotropy of the spin relaxation 
with respect to the spin polarization direction. The tem¬ 
perature and hole density dependences of the in-plane 
spin relaxation along the x-axis are plotted in Figs. [5| 
andm respectively. From Fig.jS] at the same hole density 
and temperature, the SRT with all the relevant scatter¬ 
ings included (solid curves with squares or triangles) is 
about four orders of magnitude smaller than the one with 
the intervalley hole-phonon scattering removed (dashed 
curves with dots or crosses), indicating the intervalley 
hole-phonon scattering makes a dominant contribution 
to the in-plane spin relaxation. This is very different 
from the previous results of the out-of-plane spin re¬ 
laxation, where the intervalley hole-phonon scattering is 
marginal. Specihcally, the Zeeman-like term, together 
with the intervalley hole-phonon scattering, opens an in¬ 
tervalley spin relaxation channel, as pointed out by Wang 
and Wu in ML MoS 2 )^‘^ which dominates the in-plane 
spin relaxation in p-type BL WSe 2 . 

In addition, the SRTs decrease monotonically with the 
increase of temperature (see Fig. or hole density (see 
Fig. but are insensitive to the initial spin polarization 



T(K) 


FIG. 5: (Color online) The SRTs Tsx as function of T at differ¬ 
ent hole densities. Pa = 2.5 %. Solid curves: all the relevant 
scatterings are included; Dashed curves: the intervalley hole- 
phonon scattering is removed. The inset shows the SRTs ver¬ 
sus initial spin polarization Pa at different temperatures when 

Nh = 3x10^^ cm-^ Ez = 0.006 V/A. 



(10^^ cm“^) 


FIG. 6: (Color online) The SRT Tax versus hole density Nh 
at different temperatures. Pa = 2.5 %. Ez = 0.006 V/A. 


(see the inset in Fig. El. This is due to the dominant 
intervalley spin relaxation channel, very similar to the 
case of ML MoS 2 ^ when only the intervalley electron- 
phonon scattering is considered. Specifically, the inter¬ 
valley hole-phonon scattering is in the weak-scattering 
limit, and hence the in-plane SRT Tsx = with 

TpUter representing the intervalley hole-phonon scattering 
time. Consequently, the in-plane SRT decreases with the 
enhancement of the intervalley hole-phonon scattering as 
the temperature or hole density increases. 
























B. Strong HF effective magnetic field 


1. Spin polarization dependence of out-of-plane spin 
relaxation 



FIG. 7: (Color online) The SRTs Tsz versus initial spin polar¬ 
ization Ps at (a) T = 30 K and (b) T — 300 K. Squares (Tri¬ 
angles): in the K (K') valley with all the relevant scatterings 
included; Crosses (Dots): in the K (K') valley with the inter¬ 
valley hole-phonon scattering removed. Nh = 3x10*^^ cm“^. 

E, = 0.006 V/A. 

At large spin polarization (i.e., strong HF effective 
magnetic field), different total effective magnetic fields 
[Eq. (j31)] in the two valleys show up. The out-of-plane 
SRTs in the two valleys are plotted as function of the 
initial spin polarization Pg in Figs. 0 and |8] at different 
temperatures and hole densities. We first focus on the 
case with the intervalley hole-phonon scattering removed 
(dashed curves with crosses or dots), where the out-of- 
plane spins relax independently in the two valleys. As 
seen from the figures, different SRTs in the two valleys 
are obtained, which arise from the difference of the total 
effective magnetic fields in Eq. Q. Specifically, when the 
HE and Zeeman-like effective magnetic fields share the 
same direction, say, in the K' valley, the total effective 
magnetic field 1^2^ | = |Hhf| + and hence the SRT 
(dashed curve with dots), proportional to l-f | |^ 

[see Eq. Q], increases with increasing the spin polariza- 



FIG. 8: (Color online) The SRTs Tsz versus initial spin po¬ 
larization Ps at (a) = 0.006 V/A and (b) E^ = 0.02 V/A. 

Squares (Triangles): in the K (K') valley with all the relevant 
scatterings included; Crosses (Dots): in the K (K') valley with 
the intervalley hole-phonon scattering removed. The inset in 
(b) shows the distributions for spin-down holes in the K valley 
at Nh = 2x10^® cm“^ and initial spin polarization R, = 80 % 
when t = 0 ps and t = 10 ps. T = 30 K. Nh = 2x10^^ cm~^. 


tion (i.e., HF effective magnetic field). Whereas in the K 
valley, which has opposite directions between the HF and 
Zeeman-like fields, the total effective magnetic field reads 
|r2^| = — |Hhf||, and hence the SRT (dashed 

curve with crosses), proportional to 1-|- exhibits 

a minimum at |Hhf| = \v^z \ in the spin polarization de¬ 
pendence. This minimum location in the spin polariza¬ 
tion dependence can be tuned by the temperature [see 
the comparison between Figs. [TKa) and (b)], hole den¬ 
sity [see the comparison between Fig. [TKa) and Fig. [5Ka)] 
and electric field [see the comparison between Figs. [SKa) 
and (b)]. Additionally, due to the larger effective mag¬ 
netic field in the K' valley in the present configuration, 
is always faster than . 

When the intervalley hole-phonon scattering is in¬ 
cluded (solid curves with squares or triangles), the above 
difference of the SRTs between the two valleys is sup¬ 
pressed. This is understood that the intervalley scat¬ 
tering suppresses the difference of the hole density of 
each spin between the two valleys. At low hole den- 
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sity and low temperature, the intervalley scattering is 
weak, leading to the marginal suppression on the dif¬ 
ference between and , as shown in Fig lllja) at 
T = 30 K and Nh = 3x10^^ cm“^. However, in contrast 
to this marginal suppression at T = 30 K, the suppres¬ 
sion with the same hole density at high temperature is 
markedly strong, as shown in Fig. (T^b) at T = 300 K, 
where the SRTs with the intervalley hole-phonon scat¬ 
tering included are nearly identical in the two valleys. 
This arises from the enhanced intervalley hole-phonon 
scattering. 

Moreover, even at low temperature T = 30 K, 
where the intervalley hole-phonon scattering by absorb¬ 
ing phonons is negligible, the suppression on the differ¬ 
ence between and becomes stronger with the in¬ 
crease of the hole density [see the comparison between 
Fig. [la) and Fig. [8l[a)]. This arises from the different 
spin relaxation processes in the two valleys. Specifically, 
as shown in Fig. [1] in the temporal evolution, the faster 
spin relaxation in the K valley makes the Fermi energy 
for spin-down (-up) holes larger (smaller) than the cor¬ 
responding one with the same spin in the K' valley, and 
this difference of Fermi energies [given in Eq. (IB13|) in 
Appendix IB] , 

= 0.25PsNh/Ds{e-*^^^^ - (15) 

with = Ds (-^k^Ik) ~ -^k(K') ) representing 

the density of states, triggers the intervalley scattering of 
spin-down holes from the K to the K' valley by emitting 
phonons. In this situation, the intervalley scattering time 
of spin-down holes in the degenerate limit can be written 
as [given in Eq. (ID4I) in Appendix [D] 


K(K')^K'(K) _ _ -,N 


(16) 


for > 0. Here, AA^Kl^) = with 

(C = standing for the intervalley phonon 

energy; 1/t* = 2\M^\‘^m* with \M^\ being the matrix 
element of the intervalley hole-phonon scattering. Fo¬ 
cused on a specific case at = 60 %, at low hole den¬ 
sity iV/j = 3x10^^ cm“^ shown in Fig. [71[a), the dif¬ 
ference of the Fermi energies between the two valleys 
(AP^*-'^^ |™“Ri3.3 meV) is much smaller than the inter¬ 
valley phonon energy (Hj^ j^l = 17.5 meV), and hence 
the intervalley scattering by emitting phonons is blocked 
since Nevertheless, with increasing the hole 

density to Nh = 2x10^^ cm“^, this difference over time 
(APj*-''^^|™®'’'Ril5 meV) becomes closer to the intervalley 
phonon energy, leading to the intervalley scattering by 
emitting phonons enhanced at nonzero temperature. 

In addition, we find that the above suppression on the 
difference of the SRTs at low temperature can be further 
enhanced by increasing the electric field [see the compar¬ 
ison between Fig.EDa) and (b)], which arises from a quasi 
hot-hole distribution. Specifically, as shown in the inset 
of Fig. [5Db), where the distributions for spin-down holes 


in the K valley at t = 0 ps (dash-dot curve) and t = 10 ps 
(solid curve) are plotted at the initial spin polarization 
Ps = 80 %. Due to the weak intravalley hole-phonon 
scattering at low temperature but relatively strong hole- 
hole Coulomb scattering, one finds that the fast spin pre¬ 
cessions at a large electric field (i.e., a large Rashba SOC) 
result in a quasi hot-hole Fermi distribution characterized 
by an effective hot-hole temperature Tgff (the slope of the 
dashed curve is proportional to T^) larger than T (the 
slope of the dash-dot curve is proportional to T~^), which 
enhances the intervalley hole-phonon scattering for spin- 
down holes by emitting phonons. Similar behaviors can 
also be observed for spin-up holes. With the enhanced 
intervalley hole-phonon scattering, the difference of the 
SRTs between the two valleys is further suppressed. 


2. Temperature dependence of out-of-plane spin relaxation 

The SRTs in the two valleys versus temperature T at 
different hole densities are plotted in Fig. [^ with a fixed 
large spin polarization Pg = 55 %. In the computation, 
the electric field is chosen at \rjEz\ = |Dhf| when T — 
30 K, providing a large difference of the SRTs in the 
two valleys in the absence of the intervalley hole-phonon 
scattering. 

At Nh = 0.8x10^^ cm“^ shown in Fig. [91(a), for the 
case with the intervalley hole-phonon scattering removed, 
even at the small electric field {E^ = 0.004 V/A), due 
to the large HF effective magnetic field (|Dhf|/»? = 
0.0023 V/A at T = 30 K), the total effective mag¬ 
netic field in the K' valley (|D^| = |Dhf| + hF'zl) is a 
medium field (|D^Tp| si 1) whereas that in the K valley 

(l^fffl = II^hfI - \‘nEz\\) is a small field (ID^^Tpl < 1). 
Therefore, the behaviors of the SRTs originally existed at 
small spin polarization with the small (medium) total ef¬ 
fective magnetic field [see Fig. [2Kb) ([SKb))] is observed in 
the K (K') valley here [dashed curve with crosses (dots)]. 
In addition, the difference of the SRTs in the two valleys 
becomes weaker with the increase of temperature even 
when the intervalley hole-phonon scattering is removed. 
This is because the HF effective magnetic field becomes 
weaker with increasing temperature into the nondegen¬ 
erate limit,— making the difference of the total effective 
magnetic fields between the two valleys much smaller. 
When the intervalley hole-phonon scattering is included 
(solid curves with squares or triangles), this difference 
of the SRTs is further suppressed and becomes nearly 
vanished when T > 110 K. 

At the high hole density Nh = 2x10^^ cm“^ as shown 
in Fig. [9Kb), the HF effective magnetic field (|DHF|/f? = 
0.022 V/A at T = 30 K) is insensitive to the temperature 
in the degenerate limit, and hence the SRT in the K' val¬ 
ley with the intervalley scattering removed is completely 
in regime I whereas that in the K valley sits in regime H. 
Consequently, with the decrease of Tp by increasing tem¬ 
perature, the SRT increases in the K valley (dashed curve 
with crosses) but decreases in the K' one (dashed curve 
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T(K) 


FIG. 9: (Color online) The SRT Tsz versus temperature T 
at (a) Nh = 0.8x10^^ cm“^ and Ez = 0.004 V/A and (b) 
Nh = 2x10^^ cm~^ and Ez = 0.02 V/A. The initial spin po¬ 
larization Ps = 55 %. Squares (Triangles): in the K (K') val¬ 
ley with all the relevant scatterings included; Crosses (Dots): 
in the K (K') valley with the intervalley hole-phonon scat¬ 
tering removed. The vertical black dashed and red dash-dot 
lines indicate T/ and , respectively. 


with dots). When the intervalley hole-phonon scattering 
is included, even at low temperature, the difference of 
the SRTs is markedly suppressed, thanks to the enhanced 
intervalley scattering by emitting phonons as mentioned 
above. 


3. Valley polarization 

At large spin polarization and low temperature, it has 
been understood that the faster spin relaxation in the K 
valley makes the density for spin-down (-up) holes larger 
(smaller) than that with the same spin in the K' valley, 
triggering the intervalley scattering of spin-down (-up) 
holes from the K (K') to the K' (K) valley by emitting 
phonons, as shown in Fig.[T] During this process, it is dis¬ 
covered that the initial equal occupations of holes in the 
two valleys are broken, leading to the build up of valley 
polarization. To realize the large difference of the spin re¬ 
laxation processes with , the electric field in our 



FIG. 10: (Color online) The temporal evolution of the valley 
polarization Pjv with different hole densities Nh and initial 
spin polarizations Ps at T = 50 K. 


calculation satisfies f]Ez = — DHF|t=o for given hole den¬ 
sity and initial spin polarization, and then the temporal 
evolution of the valley polarization Pn = {NK'—NK)/Nh 
at different hole densities and initial spin polarizations 
when T = 50 K are plotted in Fig. [TO] It is seen that 
over time, the valley polarization first increases and then 
decreases after reaching the maximum. This temporal 
dependence can be understood as follows. 

By using the hot-hole Fermi distribution since the hole 
distribution in the temporal evolution exhibits as a quasi 
hot-hole Fermi distribution [see the inset in Fig. |8l[b)], 
the temporal evolution of Pn in the degenerate limit can 
be written as [given in Eq. (ID5I) in Appendix [D] 


dPN 

dt 


Nh 




_K'^K V 

~p ir 


(17) 


Here, with the intervalley scat¬ 
tering time of spin-down (-up) holes is given 

in Eq. (HU), but substituting T with the effective hot-hole 
temperature of spin-down (-up) holes. 

At the beginning of the temporal evolution, 
AA^^^^|t=o = 0, and hence the intervalley scat¬ 

tering rate by emitting phonons 
[Eq. (1161) ] whereas the one by absorbing phonons is 
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FIG. 11: (Color online) The effective hot-hole temperatures 
versus the initial spin polarization Ps at different hole den¬ 
sities when t = 10 ps. Solid (Dashed) curves: for spin-down 
(-up) holes in the K valley. 


negligible at low temperature. Therefore, the out-of- 
plane spins relax independently in the two valleys with 
different SRTs and [Eq. (j!])]. Over time, due to 
the difference between and [Eq. (ITSl) ] 

increases with = AEp, and then 

becomes larger, meaning that the intervalley scattering 
by emitting phonons is triggered. Moreover, it is found 
that which arises from the larger 

effective hot-hole temperature for spin-down holes than 
that for spin-up ones. Specifically, the effective hot-hole 
temperatures for spin-down and -up holes in the K valley 
at different hole densities and initial spin polarizations 
when t = 10 ps are plotted in Fig. |TT] is found 

to increase with increasing the hole density or spin 
polarization due to the larger Fermi energy difference 
between spin-up and -down holes and the enhanced 
electric field {\riEz\ = |f^HF|t=o)- Furthermore, it is seen 
that is always larger than TJ^, which arises from 
that the spin precession at large spin polarization brings 
more spin-down holes occupying the states with the 
energies higher than the corresponding Fermi energy. 
Similar behaviors can also be observed in the K' valley. 
Consequently, with ^ ^k'^k ^ 

one has dtPN > 0 in Eq. (EZD, leading to the increase of 
the valley polarization at the first tens of picosecond. 

With further increase of time, the spin polarization 
(i.e., the HF effective magnetic field) becomes smaller 
due to the spin relaxation, leading to the decrease of the 
difference between and by reducing the differ¬ 
ence between |n^| and Then induced by 

different and at the first tens of picosecond is 
suppressed by the inter valley scattering, and hence Pn 
decreases since PnNh = Ds{AEp — AEp). Moreover, 
with Pjy > 0, AEp is relatively larger than AEp, and 


hence one has from Eq. (fT6l) since the 

difference between Tjg and Tjg (see Fig. |TT|) becomes 
smaller at the smaller spin polarization over time. Conse¬ 
quently, with the fact that the difference between 

and makes the dominant contribution in Eq. (13, 

one has dtPN < 0, leading to the decrease of the valley 
polarization after reaching the maximum. 



10 20 30 40 50 60 70 80 

Ps (%) 


FIG. 12: (Color online) The maxima of the valley polarization 
versus the initial spin polarization Ps at different hole 
densities. 


The maxima of the valley polarization P^ are plotted 
in Fig. [12] at different hole densities and initial spin po¬ 
larizations. It is discovered that this maximum increases 
with increasing the hole density and/or initial spin po¬ 
larization. This dependence can be understood from the 
following analysis. As demonstrated above, the inter¬ 
valley hole-phonon scattering suppresses the difference 
of the SRTs in the two valleys [see Fig. [TKc)]. By taking 
into account this suppression, from the KSBEs, the max¬ 
imum P^ of the valley polarization is written as [given 
in Appendix [P] 


P, 


N 


pr; ( Nj-Nj 
2 Nn \ tK 





■ ( 18 ) 


Here, the prefactor F = — e^rfff^e// 3 ^^)/ 05 ; 

tm is the time when the valley polarization reaches the 
maximum. To elucidate the trend of P^ with increas¬ 
ing initial spin polarization and hole density, we consider 
tm ~ 0 approximately since tm (around tens of picosec¬ 
ond) is much smaller than the SRT (larger than hun¬ 
dreds of picosecond). Then with = 277 and 

n^|t=o = 0 in Eq. ([4]), an estimation of P^ can be writ¬ 
ten as 


N 


F\t=t 


P T* 
^ s ‘ p 

TtK 


1 - 


1 -I- 4\r]EzTp\ 


(19) 


Therefore, with the increase of hole density or initial 
spin polarization and hence the electric field {'qEz = 
—flHF|t=o), T^z becomes faster, and P// tends to increase. 
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IV. SUMMARY 

In summary, we have investigated the hole spin re¬ 
laxation due to the Rashba SOC induced by an external 
perpendicular electric field in BL WSe 2 , with all the rele¬ 
vant scatterings included. The Rashba SOC in BL WSe 2 
is constructed from the ML Rashba SOC in each layer. 
Different from ML Rashba SOC, the out-of-plane com¬ 
ponent of the BL Rashba SOC acts as opposite Zeeman- 
like fields in the two valleys. For in-plane spins, this 
Zeeman-like field, opens an intervalley spin relaxation 
channel in the presence of the intervalley hole-phonon 
scattering, similar to the case of ML MoS 2 ^‘^ due to the 
intrinsic SOC. For out-of-plane spins, this Zeeman-like 
field is superimposed by the identical HF effective mag¬ 
netic fields in the two valleys, and hence the valley with 
same (opposite) directions between these two fields has 
a larger (smaller) total effective magnetic field. Situa¬ 
tions of weak and strong HF effective magnetic fields (i.e., 
small and large spin polarizations) are both considered. 
This difference of the total effective magnetic fields be¬ 
tween two valleys greatly enriches the physics of the spin 
and valley dynamics in BL WSe 2 . 

In the case of small out-of-plane spin polarization, the 
total effective magnetic field is determined by the 
Zeeman-like field ^MrjEz ■ We find that the intervalley hole- 
phonon scattering is marginal to the spin relaxation, in¬ 
dicating that the out-of-plane spins relax independently 
in the two valleys. Identical SRTs in the two valleys 
are obtained. We find that the behavior of the SRT is 
very sensitive to the electric field, and the different Tp 
dependences of the SRT are observed at different elec¬ 
tric fields. Specifically, the SRT can be divided into 
two regimes: regime I, the anomalous EY-like regime 
{tsz oc Tp) when IHeffTpl > 1; regime II, the normal strong 
scattering regime (tsz oc t~"^) when |Heffrp| <1. In 
regime II, the peak in the temperature (density) depen¬ 
dence due to the crossover from the degenerate to the 
nondegenerate limit is observed, similar to the studies 
in semiconductors By increasing the electric 

field (i.e., total effective magnetic field) to a medium one 
(lOeffTplRil), the crossover between regimes I and II is 
observed. At this medium field, the peak in the temper¬ 
ature dependence vanishes when the crossover between 
regimes I and II is close to the crossover from the degen¬ 
erate to the nondegenerate limit, whereas the peak in the 
hole density dependence is always observed. This arises 
from the different Tp dependences of the SRT between 
regimes I and 11. 

In the case of large out-of-plane spin polarization (i.e., 
strong HF effective magnetic field), a large difference of 
the total effective magnetic fields between two valleys is 
obtained. Due to this difference, different out-of-plane 
SRTs in the two valleys are observed when the interval¬ 
ley hole-phonon scattering is weak at low temperature 
and low hole density. Specifically, the SRT in the valley 
with the smaller total effective magnetic field is faster 
than that in the other valley. Moreover, by tuning the 


Zeeman-like and HF effective magnetic fields, the valley 
with the larger total effective magnetic field can fall into 
regime I (|HeffTp| >1) while the other valley in regime 
H (IHeffTpl < l)j s-nd hence the different behaviors of the 
SRTs between regimes I and H can be observed in the dif¬ 
ferent valleys simultaneously. This difference of the SRTs 
between two valleys can be markedly suppressed through 
enhancing the intervalley hole-phonon scattering by in¬ 
creasing temperature or hole density. Additionally, at 
low temperature in the degenerate limit, due to the weak 
intravalley hole-phonon scattering but relatively strong 
hole-hole Coulomb scattering, the fast spin precessions 
result in a quasi hot-hole Fermi distribution character¬ 
ized by an effective hot-hole temperature Tea larger than 
T, which also enhances the intervalley hole-phonon scat¬ 
tering to suppress the difference of the SRTs. 

It is interesting to discover that at large spin polar¬ 
ization and low temperature in the degenerate limit, the 
initially equal hole densities in the two valleys are bro¬ 
ken in the temporal evolution, and a valley polarization 
is built up. This arises from the different spin relax¬ 
ation processes in the two valleys at large spin polariza¬ 
tions. Specifically, in the temporal evolution, the differ¬ 
ent spin relaxation processes lead to different hole densi¬ 
ties of each spin between the K and the K' valleys, but 
with opposite signs between spin-up and -down holes. 
Therefore, the spin-conserving intervalley hole-phonon 
scattering, which transfers holes from the high density 
valley into the low density one, has opposite scattering 
directions between spin-down and -up holes. In addi¬ 
tion, the effective hot-hole temperatures for spin-up and 
-down holes are found to be different, leading to differ¬ 
ent spin-conserving intervalley scattering rates. Conse¬ 
quently, with the different intervalley scattering rates in 
the two opposite scattering directions, the initial equal 
occupations of holes in the two valleys are broken, induc¬ 
ing the valley polarization. 

Different from the out-of-plane spin relaxation where 
the intervalley hole-phonon scattering is marginal at 
small spin polarization but makes an important contribu¬ 
tion at large one, for in-plane spin relaxation, the SRT is 
found to be insensitive to the spin polarization and the in¬ 
tervalley hole-phonon scattering always makes the dom¬ 
inant contribution even at low temperature. This arises 
from the intervalley spin relaxation channel induced by 
the Zeeman-like field in the presence of the intervalley 
hole-phonon scattering, similar to the case of ML M 0 S 2 r 
and this intervalley channel dominates the in-plane hole 
spin relaxation in BL WSe 2 . 
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Appendix A: Rashba SOC in BL TMDs 

Due to the weak interlayer coupling in BL TMDs, we 
use the ML Rashba SOG^ in each layer to construct the 
BL Rashba SOG. The ML Rashba SOG by Kormanyos 
et al^ can be written as 

■f^Ra^hba = iRikxCTy - ky(J^)E^, (Al) 

where aj {j = x, y) are the Pauli matrices for the real 
spin; the ML Rashba coefficients 7 ^ can be extracted 
from the electric field dependence of the effective mass in 
the work by Zibouche et al.^ which are given in Table llll 

TABLE II: Rashba parameters of electrons 7 ^ and holes 7 ^ 
for ML M 0 S 2 , WS 2 , MoSe 2 and WSe 2 , with the unit being 



M 0 S 2 

WS 2 

MoSe 2 

WSe 2 

7h 

0.0030 

0.0326 

0.0110 

0.0416 

7r 

0.0410 

0.1072 

0.0422 

0.1166 


In BL TMDs, the Hamiltonian of the lowest four hole 
bands in the two layers near the K (K') point by Gong 
et reads 

-|- \XTzOz)l^ + tj_kTx- (A2) 

Here, At, is the spin splitting of the lowest two hole bands 
in each layer; tj_/c = —with t± and t rep¬ 

resenting the interlayer and nearest-neighbor intralayer 
hopping, respectively; stands for the Pauli matrices 
for layer pseudospin; a and A are the lattice constant 
and band gap, respectively. Th e sp ecific values of these 
parameters are taken from Ref. 

In the presence of an external perpendicular electric 
field, the total hole Hamiltonian near the K (K') point 
in BL TMDs reads 

Hy = £~\~ A^(l -t- \lTzOz )T txk'^X 

+ [r]Tz+^R{kx(Ty - kyax)]Ez, (A3) 

where the electric dipole coefficient ry can be obtained 
from the energy splitting of the lowest hole bands in the 
work by Zibouche et al.^ 

By treating the ML Rashba SOG and the interlayer 
hopping as perturbations due to the large energy split¬ 
ting Xy in each layer, we construct the effective Hamil¬ 
tonian by keeping the lowest two hole bands in two lay¬ 
ers through the Ldwdin partition methodi^i^ Up to the 


third order of the momentum, the effective hole Hamilto¬ 
nian, which is consistent with the work by Yuan et al^ 
is written as 

Hy = e^k + Ez[vk{kxcry - kyax) + mcTz], (A4) 
where Vk = v{\ -\- ak^) with v = 'y-Rt±/Xy and a = 

Appendix B: Derivation of Eq. (l4ll 

We derive Eq. (|4]) based on the KSBEs in the presence 
of the out-of-plane total effective magnetic field (along 
the z-axis) with only the long-range hole-impurity scat¬ 
tering included. In the derivation, we transform the KS¬ 
BEs into the interaction picture, and further use the 
strong-scattering approximation. 

The KSBEs in the collinear space with only the long- 
range hole-impurity scattering included can be written 
as 

^iP/ik — tRk\_kyO'x/‘^ kx<Ty /T P /2 , 

—2TTNi ^ |Vk-k'P<5(ek — £k')iPy,]i — (Bl) 

k' 

By transforming the density matrix into the interaction 
picture as 

Pyu = exp{iil'^gaz/2)pfj,]^exp{-iil'^gaz/2), (B2) 
the KSBEs in the interaction picture become 

^tPfl]^ — '^^k \_ky(Tx/2 kxCky/2j P/ik] 2'kMi 

^ ^ ( |bk—k^ ^(^k £lQ''){Pfik P^kO:(B3) 
k' 

with CTi = exp(*DggCTz/2)tTi exp(—iOgjjCTz/2). After the 
Fourier transformation 

1 

= d0kPyi,exp{-iWi,), (B4) 

one gets 

= -^fcfc/4([a+,p'-k'] - [d-JX]) - (B5) 

with 

1 m*A- 

-= —-—- d6»q|Vqp(l - cosl6>). (B6) 

Tk,l JO 

Further keeping terms |?|<1 and defining the spin vector, 
SU = Tr[^^k<^], (B7) 

one obtains 

+ (1 + \^eS^kX)dt + vlk'^rl^dt 

+ vlk'^Tk,i\SX{t) = 
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In the strong scattering limit p<Cl), Eq. (IB 8 I) 

becomes 


[2rfc,ia2 + (1 + = 0. (B9) 


With the initial condition = 0, the solution of 

is given by 


J 






t(l- v^l-c2) 
4Tfc,l/(l + |f^egTfe,i|2) 



)exp 


^(4 + \/l — cl) 

4Tfeq/(l + IflggTfeqP) 


(BIO) 


with Cz = 2y/2h'kkTk^i/{l + IflggTfeqp). In the strong 
scattering limit (c 2 <Cl), one has 

5°k.W = 5^.(0)exp(-^-^). (Bll) 

The out-of-plane SRT is therefore 

~ (1 + (B12) 

When all the relevant intravalley scatterings are included, 
by replacing with Tp, Eq. (jB12l) becomes Eq. (jH). 

Moreover, without the intervalley scattering included, 
from Eq. (IBllI) . the density difference of holes of each 
spin [Np = Tr(p^k)] between the two valleys induced by 
the different SRTs can be written into 

Ajv^( 4) ^ 0.25P,A^,,(e“‘/^- - (B13) 


Appendix C: The hole-phonon scattering matrix 
elements 

For the intravalley hole-phonon scattering, a symmetry 
group analysis of lattice vibrations at the T point, which 
belongs to the space group Deh, has been performed in 
the previous worki^ The decomposition into irreducible 
representations is as follows: 

r = Aig©2A2u©-Biu©2R2g 

©Eig©2£'iu©i?2u©2i?2g. (Cl) 

The vibration pattern of these phonon modes can be 
found in Refs. ldOl - l^ It is noted that only the modes [in¬ 
cluding in-plane OP ) and AC (Ef^) phonons 

and out-of-plane OP (4l|^, B^g) and AC (A^^) phonons], 
which induce the vibrations of the transition metal 
atoms, can trigger the hole-phonon scattering. 

Due to the weak interlayer coupling in BL TMDs^ the 
BL phonon vibrations are thought to consist of the cor¬ 
responding ML phonon vibrations in each layer, and the 
BL phonon energy spectrum is close to the correspond¬ 
ing ML one4S Therefore, for the in-plane phonons, the 
BL intravalley hole-phonon scattering matrix elements 
are constructed by using the ML ones, which have been 
reported in the work by Jin et al.^ By considering the 


number of the in-plane phonon modes in BL TMDs, the 
scattering matrix elements are given by 




MM q 




-®i>iA)g |2 


2 pULA 

{Pop? . 


(C2) 

(C3) 


But for the out-of-plane phonons, it has been shown that 
the contribution of the out-of-plane phonons in the hole- 
phonon scattering is marginal in ML TMDsi^ However, 
the out-of-plane OP phonons (A|^,R 2 g) in BL TMDs, 
which induce the relative out-of-plane vibrations of the 
transition metal atoms in the two layers, can largely influ¬ 
ence the variety of the interlayer hopping, triggering the 
intravalley hole-phonon scattering. By using the tight- 
binding model according to the arXiv version of the work 
by Viljas and Heikkila,;^ we derived the matrix elements 
of this kind of hole-phonon scattering 


|M, 


’ 2 g 12 

/i/i'ql 


P79 


2 


\M, 


42 d1 


2g |2 


t'J 2Md ^ 
2pflY' A^^ Mt ^ 


(C4) 

(C5) 


For the spin-conserving intervalley hole-phonon scat¬ 
tering in BL TMDs, with the intervalley hole-phonon 
scattering in each layer suppressed due to the large en¬ 
ergy splitting^ holes in a given valley at a given layer are 
scattered into the other valley at different layer. Only 
the inter valley phonons, which induce the relative out- 
of-plane vibrations of the transition metal atoms in the 
two layers and hence lead to the variety of the inter¬ 
layer hopping, can trigger this kind of the intervalley 
hole-phonon scattering. To derive the intervalley hole- 
phonon scattering matrix elements, one needs to know 
the intervalley phonon vibrations. As mentioned above, 
the BL phonon vibrations can be constructed by using 
the ML ones. Through the group theory analysis of ML 
lattice vibrations at the K point, which belongs to the 
space group C 3 / 1 , the decomposition into irreducible rep¬ 
resentations is as follows: 


K = 2A'©2£;[©A'©A"©A"©2A". (C 6 ) 

Only the Kq phonon modest (correspond to the represen¬ 
tation E 2 ) are the out-of-plane vibrational modes, and 
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the vibrations of the two branches of Kq phonons 
and Kq) should be two kinds of orthogonal combinations 
of the out-of-plane vibration ijjt of the transition metal 
atoms and the in-plane vibration ipd of the dichalcogenide 
atoms: 




, (C7) 

V -Mw + 

A^y2Mse/MYf■^pd - \J Mw/2Mse^t 
\/ -Mw + 2A'^Mse 


To obtain the specific combination A, one needs to solve 
the kinetic equation, which is beyond the scope of this in¬ 
vestigation. We take A = 1 approximately herei^i Then, 
we construct BL intervalley phonon vibrations by using 
the ML ones in each layer but with opposite vibration 
directions between the two layers, and by using the tight- 
binding model according to the same work by Viljas and 
Heikkila,— the scattering matrix elements are given by 



|M 


KS 


2 




(C9) 

(CIO) 


Appendix D: Valley polarization from the KSBEs 


With the hole density in each valley = Tr(p^k), by 
taking trace from both sides of the KSBEs, one obtains 


dtN^ = Tr(a*p^k|coh) + Tr(5tp^k|“r) + Tr(9tp^k|irar)- 

(Dl) 


Since the hole density in each valley is not affected by 
the spin precessions and intravalley scattering, one has 
Tr(5(p^k|coh) = 0 and Tr(i9tp^k|scl”) = 0. Additionally, 
in our calculation, the energy-splitting induced by the 
Zeeman-like term fJ.r]Ez is much smaller than the Fermi 
energy, and hence by neglecting this energy splitting in 
Eq. (O, one has 


dtN^ = 21, Y. i<,7'i'{[/;'k+q(i - /;k)(A^q+1) - /;k(i - fw+c)NM^^.+^ - 

kk'qcr/^' 

+ [/)r'k+q(l ~ /(Ik)-^q ~ /)Ik(l ~ /M'k+q)(-^q + 1)] <5(^k — ^k-l-q — (D2) 


As the hole distribution exhibits a quasi hot-hole Fermi 
distribution behavior [see the inset in Fig. |8l[b)] in the 
temporal evolution, we use the hot-hole Fermi distri¬ 
bution characterized by Tea in Eq. (ID2I) . In addition, 
one has A^q«0 at low temperature denoting 

the intervalley hole-phonon scattering through absorbing 
phonons is negligible. Therefore, with the larger density 
of spin-down (-up) holes in the K (K') valley, Eq. (ID2D 
in the degenerate limit becomes 


dN^ 

dt 


pDs 

2 


V K'^K 
P 11 


_K-s-K' )> 
> 11 


(D3) 


where 


K'(K)^K(K') ^ _-.N 

p IKH) P^ 

(D4) 


Substituting Ak' — Ak = PNNh into Eq. (ID4I) . the 
temporal evolution of the valley polarization reads 


dPN 

dt 


£l 

Nh 


^ AE^ 

( K^K' 


AAli ^ 

T It 


(D5) 


Next we derive Eq. (fT51) based on the KSBEs. Specifi¬ 
cally, it has been shown that in the degenerate limit, the 
difference of the SRTs in the two valleys is suppressed 
by the intervalley hole-phonon scattering [see Fig. EKb)]. 
With the contribution of the intravalley scatterings to 
the spin relaxation given in Eq. (IBlip . by neglecting the 
dephasing parts of the intervalley hole-phonon scatter¬ 
ing since the intervalley hole-phonon scattering is in the 
weak-scattering limit, the temporal evolution of spin po¬ 
larization in each valley becomes 


dP^ 


pp 

^ S 


D, , AE^ AE^ 






+ 


-K'^K ' 


p fr 


(D6) 


for > 0. The first (second) term on the right 

handside of Eq. (ID3p comes from the contribution of the 
spin-conserving intervalley scattering of spin-up (-down) 
holes. 


Further considering dtP^~dtPf' in Eq. (1061) . one has 


2Ds AE^ AE^ _ P^ 


jK' 


\ ^ s s 

y ~ "Zk ^k' ' 


(D7) 


‘p 4 


p 1T 
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Additionally, dtPN = 0 in Eq. (ID5I) when the valley po¬ 
larization reaches the maximum a,t t = and then by 
using Eqs. (ID4I1 and (ID7() . one obtains 


Taking and PnN^ = - APf), 

the maximum of the valley polarization [Eq. (fT51) ] is ob¬ 
tained. 


AEW) _ Pf pf 

^ ~ li57 ^ ~ ^ 


(D8) 
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